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Abstract 

Wc present the results for three-loop beta-functions for Yukawa couplings of heavy Standard Model fermions 
calculated within the unbroken phase of the model. The calculation is carried out with the help of the MINCER 
program in a general linear gauge, and the final result is independent of the gauge-fixing parameters. In 
order to calculate three-point functions, we made use of infrared rearrangement (IRR) trick. Due to the 
chiral structure of the SM a careful treatment of loops with fermions is required to perform the calculation. 
It turns out that gauge anomaly cancellation in the SM allows us to obtain the result by means of the 
semi- naive treatment of 75. 
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The Yukawa couplings being the fundamental parameters of the Standard Model (SM) Lagrangian de- 
scribe the interactions of quarks and leptons with the Higgs field. Having in mind the discovery of the 
Higgs boson candidate 0, [|[ one can hope that some day the values of Yukawa couplings will be deduced 
directly from the experimental data (see, e.g., the discussion presented in Ref. 0]). In order to obtain a very 
precise SM prediction for the running Yukawa couplings at some high energy scale, one usually uses known 
masses of quarks and leptons Q since it is the Yukawa interactions that give the fundamental fermions 
their masses after spontaneous electroweak symmetry breaking. Due to the observed hierarchy of the SM 
fcrmion masses the corresponding values arc usually defined at different scales, so one inevitably makes 
use of renormalizationgroup equations (RGE) to connect these scales. The so-called threshold (matching) 
corrections (see Refs. [f|, 0] for the case of running masses and Yukawa couplings) are also very important 
for extractions of the running SM parameters defined in the minimal subtraction (MS) scheme, in which 
counter-terms and beta-functions have a very simple polynomial structure. It also should be mentioned that 
contrary to leptons, quarks are not observed as free particles, so the pole mass which is usually associated 
with the physical mass of a particle, although being a gauge- invariant and infrared finite quantity 0, @], 
suffers from the so-called renormalon ambiguity 0, [13] • This intrinsic uncertainty of the quark pole mass is 
estimated to be of the order of Aqcd — 200 MeV. For the top-quark one usually neglects this uncertainty 
since the corresponding mass M t = 175 GcV is much bigger than Aqcd- Moreover, it is generally believed 
that due to the short lifetime the t-quark does not have time to hadronizc, so the notion of the pole mass 
can be used in this case. According to the recent studies of the relation between the MS-running mass of 
the top quark and the corresponding pole mass, electroweak corrections can become important and for the 
observed value of the Higgs boson mass can severely cancel the QCD contribution [(|. As a consequence, 
theoretical uncertainty in determination of the top-quark Yukawa coupling is reduced, thus calling for more 
precise determination of the corresponding RGEs. For all the other quarks one usually uses the MS masses 



defined initially in the context of QCD (see, e.g., Ref. [Tl| and references therein). 

The SM Higgs boson with Mh = 125 GeV decays predominately into the bb pairs. In spite of the fact 
that this decay mode is very hard to observe due to the enormous QCD background it is obvious that the 
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precise value of the corresponding Yukawa coupling is required to test whether the SM correctly describes 
Nature at the electroweak scale. If one considers leptonic decays of the Higgs boson, the most promising 
decay channel is the tree-level decay to a tau-anti-tau pair. These facts somehow motivate our study of the 
three-loop contribution to the corresponding Yukawa beta-functions. 

Moreover, we would like to stress here that a separate study of the high-energy behavior of the SM can 
also be of great importance. One can use RGE to find the scale where New Physics should enter the game, 
e.g., to unify the interactions or stabilize the Higgs potential [12 - 17 1 ■ 



One- and two-loop results for SM beta functions have been known for quite a long time [18H31| and are 
summarized in 32]. 

Not long ago full three-loop gauge coupling beta- functions were calculated [3^, HH . The beta- functions 
for the Hig gs self-coupling and top Yukawa coupling were also considered at three loops §55[. However, 
in Ref. [35|, all the electroweak couplings were neglected together with the Yukawa couplings of other SM 
fermions. In this paper, we provide the full analytical result for the three-loop beta-functions of the strongest 
Yukawa couplings yt,yb, y T for the three heaviest SM fermions (top, bottom, and tau). We take into account 
all the interactions of the SM. 

Let us briefly define our notation. The full Lagrangian of the unbroken SM which was used in this 



calculation is given in our previous paper [34| . However, we do not keep the full flavor structure of Yukawa 



interactions but use the following simple Lagrangian which describes fermion-higgs interactions 

^Yukawa = ~yt(Q^ C )t R - Vb{Q^)b R . - y T {L$)T R + h.C. (1) 

with Q = (t, and L = (v Tl r) j, being SU(2) doublets of left-handed fermions of the third generation, u R) 
t R , and t r are the corresponding right-handed counter parttQ. The Higgs doublet $ with Yw = 1 has the 
following decomposition in terms of the component fields: 

V2 



Here a charge-conjugated Higgs doublet is introduced 4> c with Yw = — 1. 
For loop calculations it is convenient to define the following quantities: 

5 ffl 92 9a Vt Vb Vr X e f t \ to\ 
3 16tt 2 ' 16tt 2 ' 16tt 2 ' 16tt 2 ' 16tt 2 ' 16tt 2 ' 16tt 2 ' ?G ' ^ G ) ' U 

where we use the SU(5) normalization of the U(l) gauge coupling g\. We also stress that the calculation is 
carried out in a general linear gauge, so the renormalization of all three gauge-fixing parameters £g> £,w, 
and £b, given in Ref. [HI, is required. 

The Yukawa beta-functions are extracted from the corresponding renormalization constants which relate 
bare couplings to the renormalizcd ones in the MS-scheme. The latter are found with the help of the 
following formulae: 



Zy f = /7 Zff 7 k 7 or Z yj . = , f = t,b,r, 



where Zffh and Zff x are the renormalization constants for the three-point vertices involving two fermions 
/ and the higgs h or the would-be goldstone boson y, respectively. The renormalization constants Zf L1 Zf R 
for left- and right-handed fermions, and Z x and Zh for the neutral components of the Higgs doublet are 
obtained from the corresponding self-energy diagrams. 

In order to extract a three-loop contribution to the considered renormalization constants, it is sufficient 
to know the two-loop results for the gauge and Yukawa couplings and the one-loop expression for the Higgs 
self-interaction. 



-'^Here we assume that neutrinos are massless. 
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The relation between the bare and renormalized parameters can be written in the following way 

oo 1 

- 2 o** = Z ak a k (ri = a k + J2ck l) ^> (5) 



Ofc,BarcM 



where pk — 1 for the gauge and Yukawa constants, pk = 2 for the scalar quartic coupling constant, and p k = 
for the gauge fixing parameters. The bare couplings are defined within the dimensionally regularized j36j 
theory with D = 4 — 2e. The four-dimensional beta-functions, denoted by /3j, are defined viafl 



ft = + ft (2) + ft (3) + • • ■ (6) 



e=0 



with pf 1 being the Z-loop contribution to the beta- function for a*. The expression for /3j can be extracted 
from the corresponding renormalization constants ([5} with the help of 

ft = E^^-^ (1) - (7) 

Here, again, aj stands for both the gauge couplings and the gauge-fixing. 

It should be noted that the divergent part of the considered three-point functions should resemble the 
tree-level structure 

AC = -^=ffh-i^=h 5 f X . (8) 

Since we separately consider the contributions to the /l/r^ and ffifL,4> vertices {4> = h, x), we have to be 
sure that the corresponding divergencies sum up to give the unit matrix in the case of the Higgs boson or 
75 in the case of %. This serves as an additional self-consistency check of our result. 

It turns out that due to the gauge symmetry the Higgs field h and the would-be goldstone boson x 
renormalize in the same way so that Z x = Z^. Moreover, the same reasoning can be applied to the 
considered Yukawa vertices giving Zffh = Zft x . This fact was also checked by explicit calculation at the 
three-loop level. 

Actually, it is not trivial to satisfy these two requirements (to conserve the chiral structure of the 
Lagrangian and do not break the gauge invariance) at three loops. Both the issues are related to the 75 
problem present in dimensionally regularized theories. It is known from the literature (see, e.g., Ref. [37| 
and recent explicit calculation |35[) that the traces with an odd number of 75 appearing for the first time in 
the three-loop diagrams require special treatment. We closely follow the semi-naive approach presented in 
Ref. (38|. First of all, we anticommute 75 with other matrices and use 7! = 1. In the "even" traces all 75 arc 
contracted with each other, so the corresponding traces can be treated naively in dimensional rcgularization. 
In "odd" traces we are left with only one 75. These traces are evaluated as in four dimensions and produce 
totally antisymmetric tensors via the relation 

tr {'f'f'fflh) = -Aie^P" (9) 

with e 0123 = 1. It is obvious that these tensors can give a nontrivial contribution only if they are contracted 
with each other. In a fully consistent calculation, this contraction should be carried out only after the 
renormalization, when one can safely take the limit e — > 0. However, one can try to use the four-dimensional 
formula 

to contract two e-tensors originating from different fermion traces before the renormalization. In (|10[) the 
square brackets denote complete antisymmetrization and 5£ is the Kronecker delta. Due to the fact that 



2 For the Yukawa couplings it is also convenient to consider the running of the coupling itself. It is obvious that /3 Vi = 
(Pi/ai) yi/2 with i = t,b, T. 
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(a) (b) (c) 

Figure 1: Sample Feynman diagrams giving rise to a non-trivial contribution due to the fermion traces 
with an odd number of 75. The second trace arises due to the projector applied to the external fermion 
line. Solid, dashed and wavy lines correspond to fermions, scalar and gauge bosons. In our calculation, we 
consider vertices with / = t,b,T and <fi = h,x- Due to gauge anomaly cancellation the contributions from 
the first two diagrams are summed up to give zero if one takes into account all the SM fermions circulating 
in the triangle loops. As a consequence, we are left with the third topology which produces only single poles 
in the regularization parameter e. 



the metric tensors in the left-hand side of Eq. (|TU|) are treated as D-dimcnsional ones, in this case the result 
of contraction can differ from the correct one only via the terms proportional to D — 4. If the considered 
diagrams produce only single poles in e, then the above-mentioned difference contributes only to the finite 
terms which we can neglect in our calculation^]. This kind of approach (semi- naive treatment of 75) was used 
in Refs. (38^ and [35| . However, if a diagram gives rise to higher poles in e, then the formal manipulations 
with antisymmetric tensors affect also lower poles and, in general, lead to a wrong result unless special 
(finite) counter-terms are introduced (see, e.g., [1^]). In order to figure out whether semi- naive treatment 
is sufficient for our calculation, we separately evaluate contributions arising from the contraction (|10[) . It 
turns out that some of the three-loop diagrams, both the fermion self-energies and vertices, give rise to the 
dangerous high-order poles in e (see Figs. fTal and ITb"|). Nevertheless, taking into account gauge anomaly 
cancellation in the SM |4(J 41 1 one can check that all these terms sum up to zero. It is worth mentioning 
that in the SM we have to put the number of colors to be equal to three (N c = 3) to satisfy the anomaly 
cancellation conditions. In the end, we are left with diagrams presented in Fig. [Tel which produce only 
single poles and give a nontrivial contribution to the final result (this kind of diagrams was also discussed 
in Ref. (H). 



As in our previous calculations we use LanHEP [421 ] and FeynArts [43j to generate relevant Feynman 
diagrams and the FORM 44[ packages MINCER 45 1 and COLOR [46| to obtain the resulting expression. Since 
MINCER is aimed only at self-energy type diagrams, we made use of the infrared rearrangement (IRR) 
trick [42| to evaluate the three-point Green functions by setting the momentum, which enters through the 
higgs external line, to zerc0- A special script was written that automatically maps the Feynarts three-point 
vertices with zero external momentum to the MINCER topologies. 

As a result of our calculation we obtain the expression^ for the three-loop Yukawa coupling beta- 
functions of the heavy SM fermions. Due to lack of space we present here only the three-loop result for the 
top quark Yukawa coupling (A = a\): 



(i) 



a t 



-8a. 



9a t 3a& 9o,2 17a% 



(11) 



3 This is due to the fact that such a complication arises only in the three- loop diagrams. 
4 A similar method was successfully used in our previous calculations l4g| . |49|| . 

5 All the results, including necessary renormalization constants, can be found online as ancillary files of the arXiv version of 
the paper. 
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(12) 



22aJa T (13) 



with n G corresponding to the number of SM generations and £ 3 = C(3)- Having in mind the anomaly 
cancellation conditions, the result is only valid for the SU(3) color group, so we substitute all the color 
invariants by the corresponding values (CV = 4/3, N c = 3, Ca = 3). 



It should be noted that the expressions are free from gauge-fixing parameters £g , £w and £,b which are 
present in the renormalization constants for the considered two- and three-point Green functions. The one- 
and two- loop corrections, given in Eqs. (fTTj) and ([T21 . arc in full agreement with Refs. 27, HI US, HH. The 
contribution (fT3)) coincides with the result of Ref. [35| in the limit of vanishing couplings 01, a 2 , a&, and a T 
(first two lines of Eq. (|T3])). 

In order to estimate the numerical impact of the calculated corrections, we use known values of the 
parameters at the Mi scale (e.g., from Ref. (HI) and compute the values of the considered three- loop 
contributions to the Yukawa coupling beta-functions. In what follows, the largest contributions that account 
for at least 99% of the total three-loop corrections to the bet a- functions are shown: 



,(3) 
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(14) 
(15) 

(16) 



where 



Oi(M z ) ' 



A 



A(M Z 



(17) 



The known results from Ref. |35| arc highlighted with boxes in Eq. (fT4l . 

It is clear that all the corrections are dominated by the terms proportional to the top Yukawa and 
the strong coupling constants. However, one can see that the a 2 a 2 contribution to the top Yukawa beta- 
function (|14j) . which was not taken into account in Ref. [35| . is comparable with the a 3 term. 

To conclude, in this paper wc present for the first time the expressions for three-loop Yukawa beta- 
functions. 



34( and [35j in a precise RGE 



The latter can be used together with the results of Refs. [3S 
analysis not only of the SM itself but also of its extensions which introduce New Physics at a scale very 
much separated from the electroweak one. However, one should keep in mind that the full three-loop beta- 
function for the quartic coupling and the anomalous dimension of the Higgs mass parameter are still missing 
and require a dedicated study which will be presented elsewhere. 
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